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Popular Summary 

Moist entropy is nearly conserved in adiabatic motion. It is redistributed rather than 
created by moist convection. Thus moist entropy and its equation, as a healthy direction, can 
be used to construct analytical and numerical models for the interaction between tropical 
convective clouds and large-scale circulations. Hence, an accurate equation of moist entropy is 
needed for the analysis and modeling of atmospheric convective clouds. 

On the basis of the consistency between the energy and the entropy equations, a complete 
equation of moist entropy is derived from the energy equation. The equation expresses 
explicitly the internal and external sources of moist entropy, including those in relation to the 
microphysics of clouds and precipitation. In addition, an accurate formula for the surface flux 
of moist entropy from the underlying surface into the air above is derived. 

Because moist entropy deals "easily" with the transition among three water phases, it will 
be used as a prognostic variable in the next generation of cloud-resolving models (e. g. a 
global cloud-resolving model) for low computational noise. Its equation that is derived in this 
paper is accurate and complete, providing a theoretical basis for using moist entropy as a 
prognostic variable in the long-term modeling of clouds and large-scale circulations. 
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Abstract 


This paper addresses an equation for moist entropy in the framework of cloud- 
resolving models. Since the energy equation is equivalent to the entropy equation with 
proper source terms, this paper for the first time derives a complete equation for moist 
entropy from the energy equation. The equation expresses the internal and external 
sources of moist entropy explicitly, providing a basis for moist entropy as a prognostic 
variable in long-term cloud modeling. In addition, an accurate formula for the surface 
flux of moist entropy from the underlying surface into the air above is derived. 

Moist entropy is used to rewrite the Neelin-Held model for the diagnosis of large-scale 
vertical velocity. After applying the model to a tropical oceanic atmosphere with mean 
annual soundings, the paper shows the sensitivity of large-scale vertical circulations to 
the radiative cooling rate and the surface flux of moist entropy, which demonstrates the 
necessity for an accurate equation for moist entropy in the analysis and modeling of 
large-scale tropical circulations. 
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1. Introduction 


The motive for an accurate equation for moist entropy has two origins: numerical 
modeling and theoretical analysis of tropical clouds and their interaction with large-scale 
circulations. On one side, the motivation comes from the rapid development of massive 
parallel computation. With the great increase in computational power, cloud-resolving 
models have significantly increased their horizontal domain size. In the near future, the 
domain size will reach thousands of kilometers so that the models can accommodate both 
clouds and large-scale circulations. With the explicit representation of clouds and large- 
scale circulations, the models can be used to explore the interaction between clouds and 
large-scale circulations. Since the models with large-scale circulations usually integrate 
for days or longer, they should represent the energy equation accurately and reduce 
computational noise, especially that from the connection between microphysical and 
dynamic processes. 

Most current cloud-resolving models use three prognostic thermodynamic variables to 
simulate non-precipitating warm clouds and those three variables plus others to simulate 
other clouds (e.g., Grabowski 1989, Tao and Simpson 1993, Tompkins and Craig 1998, 
Tao et al. 2003). The three prognostic variables are temperature (or its equivalent) and 
the mixing ratios of water vapor and cloud water. In models with the three prognostic 
variables, there are computational phenomena such as spurious supersaturation near 
cloud edges (e.g., Grabowski 1989). Special adjustment techniques are usually introduced 
to reduce their negative effects, benefiting the numerical simulation of individual cloud 
systems (e.g., Tao et al. 1989). 
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Some computational phenomena originate in the choice of prognostic thermodynamic 
variables (Grabowski 1989, Ooyama 1990, Zeng 2001). In a new trend to avoid them, 
moist entropy is used as a prognostic variable in modeling, and temperature is diagnosed 
from that and other prognostic variables (Raymond and Blyth 1986; Ooyama 1990, 2001; 
Zeng 2001). Raymond and Blyth (1986) used moist entropy and the total mixing ratio of 
airborne water (water vapor and cloud water) as prognostic variables in a parcel model. 
Ooyama (1990) analyzed the thermal relations of moist entropy and suggested moist 
entropy be used as a prognostic variable in multiple-dimensional numerical models. 
Ooyama (2001) and Zeng (2001) constructed two- and three-dimensional models with 
moist entropy as a prognostic variable, respectively, to simulate warm clouds. In their 
models neither spurious supersaturation nor negative mixing ratio of cloud water appears. 
Since moist entropy deals with water phase transition “easily”, it is logical to extend 
moist entropy as a prognostic variable to simulate cold clouds. Hence, an accurate 
equation for moist entropy is needed in the modeling and analysis of cold clouds. 

However some researchers, on the basis of the consistency between the entropy and the 
energy equations, have questioned the accuracy of moist entropy as a prognostic variable 
in long integration (Satoh 2002). In fact, the entropy equation is consistent with the 
energy equation if the sources of moist entropy (including irreversible generation) are 
expressed accurately. This paper addresses the question on the consistency between the 
energy and the entropy equations and derives accurate source terms for moist entropy 
from the energy equation. 

Likewise, there is motivation for an accurate equation for moist entropy with regard to 
research on the interaction between tropical convection and large-scale circulations. In 
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the tropics, convective forcing is thermodynamic rather than dynamic. On sufficiently 
small time- and space-scales, dynamic forcing (e.g., lifting to the level of free convection) 
is important. However, convection is self-limiting. Sustaining convective forcing depends 
upon the maintenance of favorable thermodynamic conditions. As a result, 
thermodynamics dominate on sufficiently large scales (Raymond 1995, 2000; Raymond 
and Zeng 2000). Therefore, it is imperative to represent thermodynamics accurately in 
studying the interaction between tropical convection and large-scale circulations. 

However, not all variables representing the thermodynamics approach the interaction 
straightforwardly. For example, Chamey and Eliassen (1964) used potential temperature 
to represent the energy equation to study the growth of a hurricane depression. After 
connecting the convective heating rate to the low-level convergence of water vapor, they 
introduced a concept of convective instability of the second kind. The concept is specious. 
It is questioned with moist entropy (or equivalent potential temperature) as follows. 
Assuming that the surface moist entropy is fixed and horizontally uniform, the mid- 
tropospheric temperature never surpasses the temperature of an air parcel that moves 
adiabatically from the surface to the middle troposphere, regardless of how large the low- 
level convergence is. In other words, there is an upper limit on the mid-tropospheric 
temperature. However, the upper limit is contrary to the prediction of temperature by the 
linear model of Chamey and Eliassen. 

In contrast to (potential) temperature, moist entropy is suitable for representing the 
thermodynamics involved in the interaction between tropical convection and large-scale 
circulations. Moist entropy is nearly conserved in adiabatic motion. Thus moist entropy is 
redistributed rather than created by moist convection. If a conceptual model is based on 
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the entropy equation, it can avoid implicit assumptions that address specious concepts 
such as convective instability of the second kind. Neelin and Held (1987) constructed a 
simple model of moist static energy budget, providing a starting point for understanding 
the large-scale vertical circulation in the tropics. This paper rewrites their model on the 
basis of a moist entropy equation, incorporating parcel dynamics reasoning. From the 
model, it is impossible to reach the concept of Chamey and Eliassen (1964). 

Moist entropy plays a unique role in the analysis of tropical convection. Thus, its 
accurate equation is needed in the modeling and analysis of the interaction between 
tropical clouds and large-scale circulations. However, for the time being, there is no 
accurate equation for moist entropy in a precipitating and icy atmosphere. This paper 
derives such an equation for moist entropy from the energy equation. 

The paper consists of six sections. Section 2 explains the consistency between the 
energy and the entropy equations. Section 3 derives an expression for moist entropy with 
cloud ice. Section 4 expresses the sources of moist entropy, including those in relation to 
the microphysics of clouds and precipitation. Section 5 presents a flux-form equation for 
moist entropy and derives an accurate formula for the surface flux of moist entropy from 
the underlying surface into the air above. And finally, Section 6 gives a summary and 
discusses future work. Except for special illustrations, the paper follows the symbol 
definitions in Appendix A. 

2. Consistency between the entropy and the energy equations 

This section explains the consistency between the entropy and the energy equations. 
Entropy s changes when a body is in thermal contact with its environment ( d e s ) or 
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increases as a result of internal changes within the body ( dp ). Thus, the total change in 
the entropy of a body is given by (Wallace and Hobbs 1977) 

ds = d e s + d i s (2.1) 


The change d e s is given by 



The second law of thermodynamics states 
d t s > 0 

Combining Equations (2.1) and (2.2) gives 
ds = T~ x dq + d t s 


( 2 . 2 ) 


(2.3) 


(2.4) 


Although Equation (2.3) shows the sign of dp, it can’t complete Equation (2.4) for a 
solution. The first law of thermodynamics deals with energy conservation, connecting dq 
to dT and dp. After replacing dq in Equation (2.4) with that in the first law of 
thermodynamics, the expression for djS is obtained in terms of the thermodynamic 
properties T, p and s. In other words, the entropy equation (2.4) is equivalent to the 
energy equation if d$ is represented correctly. 

On the above principle, d$ can be determined with the aid of the energy equation 
while the expression for entropy is known. Consider a dry atmosphere without diabatic 
cooling, for example, its potential temperature and entropy are 


0=niv/p)*' ,c ' 

(2.5) 

s,=C,ln(0/7; tf ) 

(2.6) 


respectively (Iribame and Godson 1981). Its energy equation is 

pdO/dt = V ■ (pK h V 6) (2.7) 
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where Kh is the eddy mixing coefficient for scalars. Dividing both sides of the preceding 
equation by #and rearranging the right side yields 

p d In 6/dt = V • (pK h V In 9) + pK h V In 6 • V In 0 
or 

pds d /dt = V • (pK h Vs d ) + pC- l K h Vs d -Vs, (2.8) 

where the last term expresses the irreversible generation of entropy due to mixing. 

In a dry atmosphere, heat passes from warm to cold areas, just as the second law of 
thermodynamics describes. The heat flux from warm to cold areas is expressed 
quantitatively by an empirical formula or the last term in Equation (2.7). It is the 
empirical formula that leads to the expression for d t s. Therefore, the empirical formula 
contains more quantitative information than the second law of thermodynamics. 

3. Moist entropy 

a. Moist entropy and equilibrium 

The concept of entropy is now clear (e.g., Morse 1969). It is discussed briefly here to 
understand different expressions for moist entropy (Betts 1973, Iribame and Godson 
1981, Raymond and Blyth 1986, Hauf and Holler 1987, Ooyama 1990, Emanuel 1994, 

Zeng 2001). A unique entropy expression makes sense only when an air parcel remains at 
an equilibrium state. However, an air parcel with three phases of water substance usually 
stays at a non-equilibrium state because saturation vapor pressure around water and ice 
particles varies with the radius and phase of the particles (Pruppacher and Klett 1978) and 
precipitating particles fall relatively with respect to air. As a result, different expressions 
for moist entropy are introduced for specific applications. 
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Moist entropy can be defined as the sum of its constituent entropies plus mixing terms 
or be introduced on the basis of entropy conservation in adiabatic motion. These two 
definitions are equivalent when an air parcel stays at an equilibrium state but are slightly 
different when a parcel stays at a non-equilibrium state. The latter definition is more 
suitable for cloud modeling and is taken in this paper. 

In fact, the issue with moist entropy in meteorology is moist entropy budget (Neelin 
and Held 1987; Emanuel 1995; Raymond 1995, 2000; Raymond and Zeng 2000) rather 
than the theoretical definition of moist entropy. Exactly speaking, the issue is the 
expressions on the internal and external sources of moist entropy while the expression for 
moist entropy is given. The next subsection derives an expression for moist entropy on 
the basis of entropy conservation in adiabatic motion. With the expression, Section 4 
derives the source terms for moist entropy from the energy equation. 

b. An expression for moist entropy 

Consider an air parcel that is in a reversible process and is thermally insulated from its 
environment. Its moist entropy is conserved. From the energy equation, variables can be 
found that are conserved. Obviously these conserved variables are either moist entropy or 
its relatives. In this way, an expression for moist entropy is obtained. 

In the approach for an expression for moist entropy, Kirchhoffs relations provide a 
starting point: 


dL,ldT = C„-c l 

(3.1a) 

dLjdT = C„-c, 

(3.1b) 

dL f jdT = c,-c j 

(3.1c) 
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where Lv, L s and Lf are the latent heats of water vaporization, ice sublimation and water 
freezing, respectively; the constants C pv , c\ and c, are the specific heats of water vapor, 
liquid and ice, respectively. The mixing ratio of water vapor q v is related to the partial 
pressure of water vapor e by 


4v = 


_R d 


K p-e 


(3.2) 


Consider a water-saturated air parcel that is in adiabatic motion and contains no 
precipitating particles. From its energy equation 


p "h ^ pvQvsw + CM< ~<lvsw ~ <li) + c i<h)dT 

- R d Td ln(p -E S J- q vsw R v Td In E sw = L f dq t - L v dq vsw 


an expression for moist entropy is derived. Using Equation (3.1) and the Clausius- 
Clapeyron equation 

dinE„/dT = L,/R,T 1 * * * (3.4) 

Equation (3.3) is rearranged to 

* + <W, ) 1” - R* 1” ^ = (3.5) 

1 ref Pref 1 1 1 ref 1 

where T ref = 273.15 K and p re / = 10 5 pa are the reference temperature and pressure, 
respectively. From Equation (3.5), the moist entropy per unit mass of dry air 1 is defined 
as 


1 Moist entropy in Equation (3.6) excludes the contribution from precipitating particles. It 

represents only the contribution from moist air that consists of dry air, water vapor, cloud water 

and ice. When cloud water and ice convert to precipitation, moist entropy passes from moist air to 

precipitating particles correspondingly. For simplicity, the transfer of moist entropy from moist 
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(3.6a) 


s = (C,+c, 9 ,)ln^-R,In-2-i+ 

‘■ref Pref 



L f (T) 

T 


<h 


where Ti is simply the air temperature for water-saturated air and is the air temperature at 
the lifting condensation level for water-unsaturated parcels that satisfies (see Section 4.b) 
(C, + C„ 9> )ln(r/j;) = (Rj + R, 9 „)to(e/£„(j;)) (3.6b) 

Using the expression for moist entropy in Equation (3.6), Equation (3.5) is rewritten 


as 


dt dt T ref T 2 dt 


(3.7) 


When neither precipitating particles nor cloud ice exists, moist entropy is conserved in 
adiabatic motion. The expression for moist entropy in Equation (3.6) is consistent with 
that of Betts (1973), Emanuel (1994) and Zeng (2001). 


4. An equation for moist entropy 

After the expression for moist entropy in Equation (3.6) is specified, the source terms 
for moist entropy are derived from the energy equation. To avoid the debate on moist 


air to precipitating particles is treated as an internal sink of moist entropy in contrast to external 
sources such as the radiative sink of moist entropy. 

The entropy of water in Equation (3.6a) is related to T ref . Thus the transfer of moist entropy 
from moist air to precipitating particles depends upon T ref . When T ref = 273.15 K, water entropy is 
small. As a result, the internal sink of moist entropy in a precipitating atmosphere is small. Its 
magnitude, as Zeng (2001) showed in his numerical model of warm clouds, is much smaller than 
that of the radiative sink for moist entropy. 
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entropy definition, s in Equation (3.6) can be defined as a “new” variable in any name 
that readers favor, and the “new” variable is used to rewrite the energy equation next. 

Let C, represent the total conversion of cloud water and ice to precipitation and E t , F t 
and S t the total evaporation, freezing/melting, sublimation/deposition of precipitating 
particles in an air parcel, respectively. The total mixing ratio of airborne water (water 
vapor, cloud water and ice) then changes as 

dqjdt = E t + S' -C, (4.1) 

Let Si represent the sublimation of cloud ice or the deposition of water vapor on cloud ice 
and Fi the freezing of cloud water or the melting of cloud ice. When air is unsaturated 
with respect to water, the mixing ratio of water vapor changes as 

dqJdt = E l +S t +S l -F l (4.2) 

where -F,- represents the melting of cloud ice to cloud water, and cloud water evaporates 
immediately in an unsaturated environment once cloud ice melts. 

a. Sources of moist entropy in a water- saturated air parcel 

This subsection derives the source terms for moist entropy from the energy equation 
when air is saturated with respect to water. A water-saturated air parcel satisfies the 
following energy equation 

(C P + C pv q vsw + c, (q, ~ q vs „ -q,+q r ) + c ; (q t +q s +q g ))dT 
- R d Td\n(p - E sw ) - q vsw R v Td In E sw (4.3) 

= L f (F t + F, - S, )dt - L s S,dt - L v dq vsw + Qdt - Q p dt 

where Q is the rate of diabatic heating and Q p the heat transfer from air to precipitating 
particles. The latter is introduced so that the precipitating particles keep the same 
temperature as the air (see Appendix C for details). Please note that -dq vsw in Equation 
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(4.3) represents the conversion of water vapor to both cloud water and ice, and -Si 
represents the deposition of water vapor on cloud ice. When cloud water exists, -S, 
measures the Bergeron process in cold clouds. 

Differentiating Equation (3.6a), then multiplying by T and combining with Equation 


(4.3) yields 


ds L f F, - LS +Q-Q T dq. ,L f qi .dr . 

— = L + c, In— -jr + - (c t q r + c, (q s + q g )))— (4.4) 

dt T T ref at T Tdt 


with the aid of Equations (3.1) and (3.4). Since the term with dT/dt is small, that term is 
calculated with an approximate value for dT/dt that is estimated from Equation (4.3). 
Neglecting the small terms in Equation (4.3) and using Equations (3.2) and (3.4) results 
in the approximation 


^ = (L, (i- +F,-S,)- L,S, + e+ Ow,„ + V)4|£) 

From the approximate hydrostatic balance, 

d\np gw 

~dt ~Fj 


(4.5) 


(4.6) 


Substituting Equations (4.1), (4.5) and (4.6) into (4.4) gives 

~ = — ^ + ' Q —- + c i(S t -C t )\n-^— + (L f q i -T(. Cl q r +c t (q s +q g ))) 


dt 

f 


l A v . 


re f 

V 


Lf(F i +F l -S i )-L s S l +Q-(l+^)gw 




R d T 


J 


(C p + C pv q )T i + 


2 tfq ^ ' 

l i vrlvsw 


(4.7) 


K 


J 


where the right side represents the source of moist entropy when air is saturated with 
respect to water. 
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b. Sources of moist entropy in a water-unsaturated air parcel 

When an air parcel is unsaturated with respect to water, its energy equation is 

(C p + C pv q v + c,q r + c, (q t +q s +q g ))dT- R d Td In (p -e)- q v RJd In <? 

= L f F,dt-L v E l dt-L s (S,+S i -F i )dt + Qdt-Q p dt ' } 

where -F, represents the melting of cloud ice to cloud water, and cloud water evaporates 
immediately once cloud ice melts. 

Before deriving the source terms of moist entropy from the preceding equation, some 
thermal relations are needed for a water-unsaturated parcel. First is an expression for 7/ 
the temperature at the lifting condensation level. If a water-unsaturated air parcel is lifted 
adiabatically to the lifting condensation level, water vapor and energy are conserved. 
That is 

to(«/£„(T, ))=!«((/> -«)/(p, -£„(r,») (4.9a) 

and 

(c, + C pv q t )ln(T/T l ) = R d 1 4(p-e)/(p, - £„( 7))))+ K A ln(e/£„(7;)) (4.9b) 

where pi is the total air pressure at the lifting condensation level. Eliminating pi in 

Equation (4.9) yields one equation for 7) or Equation (3.6b). 

From Equation (3.2), e is represented with p and q v as 

dine _ dhi p | R d _ dq, (4 1Q) 

dt dt q v (R d +R v q v ) dt 

Differentiating Equation (3.6b) and substituting dT into Equation (4.8) yields the 
tendency of the temperature at the lifting condensation level 
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din 7] LjF'-LyE -L (S, + Sj-F^+Q-Q dT 


dt 


Tdt 


_( + (/j + ^ _(^ )i n Z.)^r)((7 +(7 ^ + R v c lv A>(^t) \-i 

(4.11a) 

with the use of Equations (3.4) and (4.10), where the symbol q is introduced as 

q=q v + (q, -q v ~q i +q r )+ c, (q, +q s +q g ))/ c„ (4. l lb) 

The source terms for moist entropy are derived from Equations (3.6) and (4.8). 
Differentiating Equation (3.6a), then multiplying by T and combining with Equation (4.8) 
yields 

ds L f F l —L v (E t + S'- F^-L'S'+Q-Q T dq t 

— = — — + c, In — 

r r > i rr> j. 


dt 


T re j dt 


+ q v (Cpy ~ c i~ + (4, + (c f - c z )?, - c P vq)j^ t ( 4 - 12 ) 

t Z v (7j) i dg v dlnp 

7) R d +R v q v ’ dt dr 

with the aid of Equations (3.4) and (4.10). Since the terms with dTIdt and dkip/dt are 
small on the right sides of Equations (4.11a) and (4.12), they are approximated using 
Equation (4.6) and 


dr rj din p + LfIiZ hEi ~ h^i + 5 - zIjIIR 


dt C p + C^q dt 


C p + C pv q 


(4.13) 


which is obtained from Equation (4.8) after the small terms are ignored. Combining 
Equations (4.1), (4.2), (4.6), (4.11), (4.12) and (4.13), the source of moist entropy is 
expressed as 
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ds_ _ L f F t - L s S, +Q-Q i 
dt 


~ * ' — +Cf(E t +S t -C,)ln-f— 

l rej 


+ <7v (C pv - c, - + (4, + c,<7, + (c, - c, )g, - C pv g*)y 


(4.14) 


+ + — ^—)s, + (A^ZI) .A. + + s t - F t ;.) -^-q v w 

T t R d +R,q v ‘ T t T R d +R v q v A • '' R d T H 

where the new symbols 


LjFt-Lfet-mS'+St-FJ + Q-gw 

* C p +C pv q 

^ _ ( I,F, -LA -ms, +S t , -F,) + Q-C„(g*-g,H- 


(4.15a) 


(■ 


K 


(R d +R v q v )q v R d +Kq v 

(C +C n + 


->•(«. +s, +«,-/;)) 


7\ 


(4.15b) 


The approximations in Equations (4.5), (4.6) and (4.13) introduce very small errors in 
the moist entropy source terms. Scale analysis (omitted) shows that the errors are much 
smaller than the moist entropy sink from the radiative cooling of l°C/day. 


5. Moist entropy budget 

a. Flux form of the moist entropy equation 

The equation for moist entropy in a precipitating and icy atmosphere is written in a 
flux form as 

3,/* + V -(pvs) = p(R s +M S +G S +D S ) (5.1) 

where R s , M s , G s and D s denote the sources of moist entropy from radiation, microphysics 
of clouds and precipitation, generation due to mixing, and turbulent diffusion, 
respectively. The moist entropy radiation term R s is defined as 
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R S =Q/T 


(5.2) 


From Equations (4.7) and (4.14), the microphysics term M s is expressed as 


M 


= — t' S ' Q, +c,(S, -C,)1n^-+(L,q, -T(c,q, + <:,(«,+«,») 


l ref 


■^V^V 


v 


L f (F i +F l -S i )-L s S l +Q-(l+^)gw 


R d T 


when air is saturated with respect to water, and 

M, = L,F, ~ L f' ~ Q? +c,(E, +S, -C,)ln-L 

1 l ref 


(C,+C^q')T 2 + 


L v q v 


, (5.3a) 


R 


V J 


+ <7v (C pv - C, - + (9, + c,$ ( + (c, - c, )$ - c pv q*)^- 

+ ^ ) 5 , + (Mil— )(£ + 5 w 


(5.3b) 


when air is unsaturated with respect to water. For simplicity, M s represents a small part of 
the radiation contribution through Q and Q in An- 

Turbulence or mixing brings about the diffusion and generation of moist entropy. The 
diffusion term D s is expressed as 


D, = p v • (pKJJs) 

From Equations (2.8) and (B2), the generation term due to mixing G s takes the form 


(5.4a) 


G 5 =c; l K,ys d -Vs d 


K h , 

+ y-(c,-c 1 - T 


>L(~ Lf (T) ^^ V7T 1 

V-pv 


)v 9 , ■ vr -&-(c„ - c, - iSi)v 9 , • vt; 

/ 1 1 


(5.4b) 


+ p“‘V-p^ 


f (c, - c, -^)vr -^(c„ -c, -^)vr, 


Mi 

r 


& 

5 


Mi). 

2; 


ft. Surface flux of moist entropy 
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The surface flux of moist entropy from the underlying surface into the air above plays 
an important role in the moist entropy budget. This subsection derives its expression in 
terms of Qh and Qe the sensible and the latent heat fluxes from the underlying surface to 
the air. 

Consider a thin layer with depth h adjacent to the underlying surface. When air in the 
layer is saturated with respect to water, the sensible and the latent heat fluxes increase 
temperature and water content as 


(r , r p-dta (P-EJ 3 ln£ JW Q H +Q E 

{C ' +C > > ir R ' r — s — 


(5.5a) 


dq, _ Qe 

dt phLJT) 


(5.5b) 


Since the air is saturated, the mixing ratio of water vapor changes with temperature as 


MT) n ar 

dt qvswK R V T J Tdt 


(5.6) 


which is derived from Equations (3.2) and (3.4). Differentiating Equation (3.6a) with 
respect to T and multiplying the resulting equation by ph yields the increase of moist 
entropy in the layer or the surface flux of moist entropy from the underlying surface to 
the air 


F s = Cl 


T^trln j-*iC r + C[ q, -R d +q vlw (C -c, +(M 2 - 2 )M^) 

A- CO T ref R v T 1 

_ c 

*' T T(C p + C pv q -R d - q vsw T~ l L v (D) 


(5.7) 


with the use of Equations (5.5) and (5.6). 

When air in the layer is unsaturated with respect to water, 
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(5.8a) 


(C +r Y)— fl pr 9lne 

( C p +c PV q ) dt R d dt 9v J dt 



fyy = Qe 
dt phL v (T ) 


(5.8b) 


Differentiating Equation (3.6b) with respect to T and making use of Equations (3.2), (3.4) 
and (5.8), the tendency of temperature at the lifting condensation level is obtained as 


f 
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Differentiating Equation (3.6a) with respect to T and then multiplying by ph yields the 
surface flux of moist entropy 
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(5.10) 


with the use of Equations (5.8) and (5.9). 

Equations (5.7) and (5.10) express accurately the surface flux of moist entropy from 
the underlying surface into the air. They are used, as an example, to calculate the surface 
flux of moist entropy over ocean. The sensible and the latent heat fluxes over ocean take 
the following bulk parameterization schemes 
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Q H =PC p c d W s (T su -T) 

Q E =pL v c <l w s (q YSS -q v ) 


(5.11a) 


(5.11b) 

where the drag coefficient c^l.lxlO' 3 ; W s is the surface wind; and q vss is the saturation 
mixing ratio at T sst the sea surface temperature. 

The surface flux of moist entropy is usually estimated, for simplicity, with the 
following approximate formula 

F sappr = P C d W s ( S sea ~ S air ) (5-12) 

where s sea and s air represent the moist entropy of the surface sea and air, respectively. This 
approximate formula can be compared with the accurate formula in Equation (5.10) to 
show its error. Figure 1 displays the error quantities F s -F sappr and ( F s -F sappr )/F s versus 
surface relative humidity, when T ss f= 302 K, surface pressure p=1016 hpa, surface air 
temperature T= 299 K, surface wind Ws=5 m s' 1 , and there is neither ice nor precipitating 
particles. Figure 2 is the same as Figure 1 except that the air surface moist entropy 
s ai r= 230 J kg' 1 K' 1 instead of T is fixed. As both figures show, the surface flux of moist 
entropy in Equation (5.12) contains about 2% relative error that is equivalent to a moist 
entropy sink from atmospheric cooling of 0.04°C/day. Such error, as shown in the next 
subsection, can distort the modeling of large-scale vertical circulations in the tropics. 

c. The moist entropy budget in a tropical atmosphere 

Using the equation for moist entropy, this subsection analyzes the moist entropy 
budget in a tropical oceanic atmosphere to show the necessity for an accurate surface 
moist entropy flux in modeling large-scale vertical circulations. Neelin and Held (1987) 
constructed a simple model on the moist static energy budget, providing a starting point 
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for understanding large-scale vertical circulations in the tropics. In contrast to moist static 
energy, moist entropy is a conserved variable for an air parcel in adiabatic motion. Thus, 
their model is rewritten in moist entropy so as to incorporate a parcel dynamics 
methodology. 

The internal source of moist entropy is very small in contrast to F r , the radiative flux 
of moist entropy at the tropopause, and F s , the surface flux of moist entropy from the 
ocean into the air above (Zeng 2001). Thus, in a quasi-steady atmosphere, the inflow of 
moist entropy is balanced by the outflow of moist entropy, or 

F,~F r = f -\ h g' l dp (5.13) 

J p, 

where the horizontal difference in moist entropy y h -Vs is ignored (Neelin and Held 1987). 
Just as Equation (5.13) states and Figure 3 shows schematically, the increase in moist 
entropy through the bottom of the atmosphere and the top must be balanced through the 
lateral boundary by the large-scale vertical circulation. Hence the large-scale vertical 
circulation is controlled by F s -F r with the aid of the mass continuity equation 

f‘V^> = 0 (5.14) 

J Pt 

The model in Equations (5.13) and (5.14) is applied to the atmosphere with the 
average tropical profile (Jordan 1958) to diagnose large-scale vertical velocity. The 
average profile is displayed on the left side of Figure 3. Thick and thin solid lines 
represent moist entropy and saturation moist entropy versus height, respectively. Points A 
and B indicate surface moist entropy and saturation moist entropy, respectively. Thick 
and thin dashed lines represent the moist entropy and saturation moist entropy of a lifted 
surface air parcel, respectively. When the surface air parcel moves upward, it passes the 
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lifting condensation level (LCL) first, then the level of free convection (LFC), and finally 
reaches Point C with zero buoyancy. Point C usually corresponds to the tropopause. Let 
p t represent the pressure at Point C. Thus large-scale vertical circulations are confined 
below the height p t . 

Figure 3 shows that moist entropy decreases with height at first, and then increases 
with height. The vertical profile of moist entropy is fitted with the following function 




P s P s mil 

. Pt P s min 


(P-Psmm ) 


Psmst < P < Ps 
P,<P< Psm in 


(5.15) 


where p s is surface pressure, and p S mm~7 00 hpa is the pressure with the minimum moist 
entropy s m i„. In the atmosphere with the preceding vertical distribution of moist entropy, 
small cumulus clouds and large turbulent moist eddies lead to the down-gradient 
transport of moist entropy between p s and p sm m, and deep convective clouds lead to the 
up-gradient transport of moist entropy between p smin and p, (Riehl and Malkus 1958). 
With buoyancy sorting (Raymond and Blyth 1986), cumulus clouds lead to a large-scale 
vertical velocity with a maximum value at the level p wmax higher than p sm m, which 
corresponds to the concept of “gross moist stability” that Neelin and Held (1987) 
proposed. Thus it is assumed that the vertical distribution of horizontal divergence 


V-v A = 


P Pw max 
Ps Pw max 
P_ Pwmax 
S-Pt Pwmax) 


2 (-Ps Pw max)^ ' ^ hs 


Pwmax < P < Ps 
Pt < P < Pw max 


(5.16) 


where V-v^ is surface horizontal divergence and pwma*=650 hpa is estimated. 


21 


Equations (5.15) and (5.16) determine the net lateral output of moist entropy that is 
caused by the large-scale vertical circulation. Substituting Equations (5.15) and (5.16) 
into (5.13) yields a solution of the surface horizontal divergence V-v fc . Then the vertical 
distribution of large-scale vertical velocity is obtained as a function of F s -F r , s min , p smi „ 

and Pwmax- 

In the model in Equations (5.15) and (5.16), the large-scale vertical circulation is 
controlled by F s -F r . Since the tropical atmosphere stays close to radiative-convective 
equilibrium, F s is close to F r . Therefore (F s -F r )/F s is very small. If the formula for the 
surface flux of moist entropy such as Equation (5.12) contains ~2% relative error, it 
distorts the modeling of the large-scale vertical circulation. Summarily, Equations (5.7) 
and (5.10) express the surface flux of moist entropy accurately. They are therefore 
needed for long-term cloud modeling, where large-scale vertical circulations are involved. 

6. Summary and future work 

This theoretical paper deals with the equation for moist entropy in a precipitating and 
icy atmosphere. On the basis of the consistency between the energy and the entropy 
equations, the paper derives accurate source terms for moist entropy from the energy 
equation, including those in relation to the microphysics of clouds and precipitation. As a 
result, the paper obtains a complete equation for moist entropy with explicit internal and 
external sources. Since the equation is equivalent to energy equation, it provides a basis 
for using moist entropy as a prognostic variable in the next generation of cloud-resolving 
models (e.g., a global cloud-resolving model). 
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In addition, the paper derives an accurate formula for the surface flux of most entropy 
from the underlying surface into the air above. Such an accurate formula is necessary in 
long-term cloud modeling, because the tropical large-scale vertical velocity is sensitive to 
the surface flux of moist entropy. 

Moist entropy is nearly conserved in adiabatic motion. It is redistributed rather than 
created by moist convection. Thus an accurate equation for moist entropy, as a healthy 
direction, can be used to construct analytical and numerical models of moist convection. 
In this approach for the next generation of cloud-resolving models, moist entropy will be 
used as a prognostic variable, and this paper provides a theoretical basis. In addition to 
this work on entropy budget, other numerical techniques need to be developed to treat 
sound waves and diagnose the temperature for correct energy and entropy conversion. 
Those numerical techniques will be introduced elsewhere. 

Analytical models of moist convection written in moist entropy, just as that in Section 
5.c, are usually based on the assumption that the internal sink of moist entropy is small 
enough to be ignored. Zeng (2001) used a three-dimensional numerical model of warm 
clouds to show that the internal sink of moist entropy is very small in contrast to the 
moist entropy sink from the atmospheric cooling of l°C/day. Since ice plays a role in 
some atmospheric phenomena, the assumption needs to be re-checked when ice is 
involved. This paper provides an expression for the internal sink of moist entropy. The 
expression will be used to estimate the magnitude of the internal sink of moist entropy 
with the aid of cloud-resolving modeling. 
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APPENDIX A 


List of Symbols 

Cp/Cpv/c/Ci : specific heat of dry air/water vapor/liquid water/ice 
C, : total conversion rate from cloud water and ice to precipitation 
D s : subgrid diffusion of moist entropy due to mixing, see (5.4a) 
e : partial pressure of water vapor 
E m : saturation vapor pressure over water 

E/F/St : total rate of evaporation/freezing/sublimation of precipitating particles 
Fi : rate of cloud water freezing or cloud ice melting 
F r : radiative flux of moist entropy at the tropopause 

F s : surface flux of moist entropy from the underlying surface into the air above, see (5.7) 
and (5.10) 

g : acceleration due to gravity 

G s : generation of moist entropy due to mixing, see (5.4b) 

Kh : eddy mixing coefficient for scalars 

EJLs/Lf'. latent heat of vaporization/sublimation/freezing, see (3.1) 

M s : microphysical source of moist entropy, see (5.3) 

p : total pressure of moist air 

p re f= 10 5 pa : reference pressure 

q : heat absorbed by a unit mass of air 

q * : defined in (4.11b) 

q/qc/q/q/q/q g : mixing ratio of water vapor/cloud water/cloud 
ice/rainwater/ snow/graupel 
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q vsw : saturation mixing ratio of water vapor over water 
qf=q v +q c +qi : total mixing ratio of airborne water 
Q : rate of diabatic heating 

Qh/Qe '■ sensible/latent heat flux from the underlying surface into the air above 
Q p : heat transfer from air to precipitating particles, see (Cl) 

Rd/Rv : gas constant of dry air/water vapor 
R s : radiative source of moist entropy, see (5.2) 
s : (moist) entropy of a unit mass, defined in (3.6) 

Sd : entropy of dry air, see (2.6) 

Si : sublimation of cloud ice or deposition of water vapor on cloud ice 
t : time 

T : temperature 

Ti : air temperature at the lifting condensation level, defined in (3.6b) 

T re f =273. 15 K : reference temperature 

v : velocity vector 

Yh : horizontal velocity 

VrJVJV g : fall speeds of rain/snow/graupel 

w : velocity velocity 

z ' height 

p : air density 

Aj/A-n : defined in (4.15) 
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APPENDIX B 


Irreversible Generation of Moist Entropy due to Mixing 

Section 2 derives the irreversible generation of dry entropy due to mixing. This 
appendix discusses the irreversible generation of moist entropy due to mixing. For a 
moist atmosphere, there are the following two equations 


p ^ + ... = V . (pKhWqv) 


(Bla) 


p-£-+- = V-(pK h V qi ) 


(Bib) 


where the right sides represent the effects of mixing. Multiplying Equation (Bla) by 
Lv(Ti)/Ti, subtracting (Bib) times Lj(T)/T, and rearranging terms yields 

p h2±*k- p hv± A + ....v.(*,v(£Si, -Ml,.), 

T, dt T dt h Tl Hv T q,)) 

-c, ■ vr-^*-(c_ -C, -MUyvq^T, (B2) 


+v -P K l^r(fi -C, 

K 1 1 h h 


with the use of Equation (3.1). On the right side of the above equation, the last two terms 
represent the irreversible generation of moist entropy due to mixing except for the last 
term in Equation (2.8). 
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APPENDIX C 


Energy Transfer from Air to Precipitating Particles 

Precipitating particles fall with respect to air. As a result, they move downward from 
one parcel to another parcel. Since their temperature is usually lower than the 
atmospheric temperature, heat is transferred from the air to the precipitating particles. 
This appendix derives an expression for the heat transfer from air to precipitating 
particles. 

In a Cartesian coordinate system moving with the air, a cell with unit horizontal area 
spans the height z to z+Az. Precipitating particles fall into the cell at the top with the flux 
p(z + kz)q j(z + Az)Vj , where V} denotes the terminal velocity of precipitating particles. 
Assume that the particles distribute uniformly in the cell and keep the same temperature 
as the air. Thus, the particles change their temperature from T(z+Az) to T(z). As a result, 
the precipitating particles absorb energy from the air at the rate 

1 dT 

CjPqjVj (T(z + Az / 2) - T{z + Az))/Az or — CjPqjVj — — in the unit volume, where cj is 

2 oz 

the specific heat of either water or ice. 

Precipitating particles are rainwater, snow and graupel/hail. Thus the heat transfer 
from air to precipitating particles in unit mass of dry air is expressed as 

Q, =-\{WrV,+c,q,V,+ c ,<l s V,)^ (Cl) 

In other words, precipitating particles bring about energy exchange between air parcels. 
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Figure Captions 


Figure 1 Error in the surface flux of moist entropy in Equation (5.12) with respect to 
surface relative humidity when the air surface temperature T-299 K is fixed. 


Figure 2 the same as Figure 1 except that the air surface moist entropy s air =230 J kg' 1 K 
1 is fixed. 


Figure 3 Vertical distribution of moist entropy (thick solid line) and saturation moist 
entropy (thin solid line) and a schematic of the large-scale vertical circulation in the 
Tropics. Points A and B indicate surface moist entropy and saturation moist entropy, 
respectively. Thick and thin dashed lines represent moist entropy (Line A-LFC-C) and 
saturation moist entropy (Line B-LCL-C) of a lifted surface parcel, respectively. 
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Error of Surface Entropy Flux (J/K/s/m 2 ) 



Figure 1 Error in the surface flux of moist entropy in Equation (5.12) with respect to 
surface relative humidity when the air surface temperature T = 299 K is fixed. 
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Radiation Flux of Moist Entropy 



Figure 3 Vertical distribution of moist entropy (thick solid line) and saturation moist 
entropy (thin solid line) and a schematic of the large-scale vertical circulation in the 
Tropics. Points A and B indicate surface moist entropy and saturation moist entropy, 
respectively. Thick and thin dashed lines represent moist entropy (Line A-LFC-C) and 
saturation moist entropy (Line B-LCL-C) of a lifted surface parcel, respectively. 
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